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Abstract. A rudimentary calculation is employed to evaluate the possible effects 
of /3-decays of excited state nuclei on the astrophysical r-process. Single particle 
levels calculated with the FRDM are adapted to the calculation of /3-decay rates of 
these excited state nuclei. Quantum numbers are determined based on proximity to 
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Nilson model levels. The resulting rates are used in an r-process network calculation 
in which a supernova hot-bubble model is coupled to an extensive network calculation 
including all nuclei between the valley of stability and the neutron drip line and with 
masses 1<A<283. /3-decay rates are included as functional forms of the environmental 
temperature. While the decay rate model used is simple and phenomenological, it is 
consistent across all 3700 nuclei involved in the r-process network calculation. This 
represents an approximate first estimate to gauge the possible effects of excited-state 
/3-decays on r-process freezeout abundances. 

PACS numbers: 26.30.Hj,26.50.+x 
1. Introduction 

The r-process is responsible for the synthesis of roughly half of all nuclei heavier than 
A~70 and all of the actinides [U [2] . The solar system r-process abundances act as the 
canonical constraint to r-process theories as well as the prime indicator of the success 
of r-process models. Several r-process sites have been proposed; the hot-bubble region 
of a type II supernova (SNII) has been modeled fairly successfully. The composition of 
the environment in which the r-process occurs might be expected to have a profound 
effect on the final abundance distribution. Observations indicate that the r-process site 
is primary [3] and further evidence may suggest that the r-process is also unique [3]; 
it may occur in a single site or event. The uniqueness of the r-process site, however, 
remains a subject of study [5]. 

Nuclear properties also constrain the r-process, and the purpose of this work is the 
examination of one particular characteristic - /3-decay - as it relates to the r-process. 
The /5-decay inputs, and other nuclear physics inputs, have been shown [6j [7] to have 
important effects on the success or failure of r-process models. This is somewhat 
unfortunate, as properties of only a few nuclei on the neutron closed shells closer to 
stability have been experimentally determined, while data for the rest are relegated to 
calculation pp. Of paramount importance is the determination of nuclear masses and 
/3-decay rates. Nuclear mass formulae based on the microscopic properties of nuclei are 
slowly replacing the empirical droplet models, and these can change resulting reaction 
rates by factors as large as 10 8 [8]. As well, the r-process path is affected by the choice 
of mass formula, since the path roughly follows a line of constant S n [1] . 

A successful r-process calculation must predict an abundance peak for nuclei in the 
A~195 region - a difficulty over much of the r-process parameter space. However, as 
discussed below, /3-decays from excited state nuclei may help to mitigate this difficulty. 
They do so by allowing the r-process to proceed at a faster rate, thereby enhancing the 
abundances at higher masses. 

For the purposes of this study, the most recent semi-gross theory of /9-decay [9] has 
been adapted to neutron-rich nuclei relevant to the r-process. The ability of this model 
to determine decay properties of an extremely wide range of nuclei with reasonable 
accuracy and speed makes it ideal for this preliminary calculation. In particular, the 
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semi-gross theory has good agreement for very neutron-rich nuclei [TUl [TT] . It has 
also been used to improve the accuracy of decay rates for astrophysical calculations by 
incorporating first-forbidden transition strengths [12]. In its original form, the gross 
theory of /5-decay assumed that the energy states of a nucleus consist of a smoothed 
distribution with transition strengths that peak at or near the energy of the isobaric 
analog state [13], [H] . Subsequent evolutions of the gross theory incorporated strength 
functions allowing for transitions of higher forbiddenness [15] , as well as improvements 
over the original theory to include odd-odd effects [16], sum rules [17] . even-odd mass 
differences |18] . and improvements on the strength functions |19j . 

While the precision of the Fermi gas model used is consistent with that of the gross 
theory itself, more accurate models might be used. Recently, shell effects in the parent 
nucleons have been taken into account by using an energy distribution for single-particle 
states [9l ED]; this is denoted as the "semi-gross theory." In these models, the energy 
distribution of the daughter states is still assumed to be smooth. Thus, the transition 
strength functions depend on the quantum numbers only of the initial parent states 
and are independent of the energy of the parent state; transition types are based on a 
statistical weight for a particular parent state to make a particular type of transition. 
The advantage of using single-particle strength functions is that quantum numbers can 
be assigned to the states easily, lending a better notion of the actual strength of each 
type of transition involved. 

Since decay rates of nearly all of the nuclei along the r-process path have yet to 
be studied in a laboratory, r-process calculations rely heavily on calculated decay rates. 
Further, the temperature of the r-process environment ( 10 9 K) necessitates accounting 
for nuclei in excited states, especially given the expected high level density of these 
far-from-st ability nuclei. Some of the effects that might be expected if one considers 
excited state nuclei in r-process simulations include increased (n,7) and (7,n) rates, 
which might shift the r-process path, but would tend to counteract each other as the r- 
process is generally presumed to proceed at (n,7)<->(7,n) equilibrium, increased neutrino 
spallation rates, tending to enhance smoothing in post-processing, and increased beta- 
decay rates. However, at signicant excitation, neutron separation energies are low 
enough that neutron emission may be a dominant decay mode. This is discussed briefly 
along with the discussion of excited-state decays in the network calculation. 

Section El of this paper is an overview of /3-decay calculations and average properties 
in calculating decay rates calculated using the more recent semi-gross theory. This 
includes a brief review of the semi-gross theory of /?-decay and a description of the 
calculations of single-particle states and their relationship to /3-decay in §2.11 Energy 
levels are calculated in §2.21 using the Finite Range Droplet Model (FRDM). Results of 
/3-decay calculations of excitated state nuclei with the adapted FRDM are discussed in 
§2.31 along with their potential astrophysical importance. Application of these results 
to an r-process network calculation is made in §|3] with results from a model with 
several environmental parameter sets described in £JU Future work and possibilities 
are discussed in £j5j along with experimental possibilities. 
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2. Average Nuclear Properties in /3-Decay 

Average theories of /3-decay, such as the gross theory and the more recent semi-gross 
theory, make them useful in estimating and adapting decay rates of several nuclei to 
large systems such as astrophysical network calculations. The major assumption in 
these theories is that nuclear level densities are high enough to make integration over 
states a fairly accurate approximation by replacing the matrix elements with transition 
probability distribution functions. For example, the decay rate for Fermi and GT 
transitions, as calculated in the gross theory by Takahashi and Yamada, is given by 
P-D: 



A 



m 5 c 4 r o rei 



2nh 7 



/ C £l dN 
/ Y j D n {E,e)fa{-E)-^W{E,s)dedE (1) 
QJe {E) n « e 

which integrates over the product of the transition probability function Dq(E,s) 
between a parent nucleus state at energy e and a daughter state with a transition 
energy E, the density of states dN\/de of the decaying nuclei; the weighting function of 
the final state nucleus W(E,e), which imposes the Pauli exclusion principle on the rate 
calculation; and the form factor fn(-E), which is the product of the emitted electron 
wave functions and the Fermi function for a transition of type Q. Terms used in the 
equations of this paper are summarized in Tables [TJ - [2j 



The density of parent states supplies information on the nucleons available for 
decay, while the weighting function controls the availability of the daughter states to 
which the nucleons can decay. In the case of an even number of daughter nucleons, the 
available energy levels are filled up to a value £i-Q, and in the case of an odd number, 
there is a hole at the highest energy in the ground state. Therefore, the weighting 
function is given by [T4] : 



W(E,e) 



e + E>e!-Q 
e + E < ex - Q 

(2) 

1 e + E>e 1 -Q + 2A 

U, { £ e + E = e 1 -Q 

otherwise 

The value of n3, the number of daughter nucleons promoted to the highest energy due 

to the pairing forces, is: 

re 1 ~Q+A , d j^ 

de 

where dN 2 /de is the density of states of protons (neutrons) in (5~^ decay and A is the 
pairing gap. 

The decay must be energetically possible, so the lowest possible energy that a 
decaying nucleon can have is £i-Q-E, recalling that E is the (negative) transition energy 
of the decaying nucleon. Thus, e in Equation [1] equals £i-Q-E, removing that variable 
from the equation. 



n 3 = I i^)de (3) 
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2.1. Single- Particle States 

Although a smooth model is adaptable to decays, particularly from the ground-state, 
excited states present the added difficulty of potentially altering the strength functions. 
For an accurate treatment, an estimate of the quantum numbers of these states is made. 
In the initial calculation, discrete energy levels are used. Thus, the density of states for 
neutrons (protons) in [3~^ becomes: 

(£) w/ , = 5>(<-«.> w 

where the value of n, corresponds to the number of nucleons in each state i. For discrete 
Fermi states, the levels are two-particle levels with n=zero, one, or two. The weighting 
function W(E,e) in Equation [1] can also be adjusted for single particle daughter levels. 
In the original gross theory, the weighting function takes on values between zero and 
one and is typically zero or one except at a few special points. If energy levels are 
discrete, then the value of the weighting function for a two particle level is either zero, 
one, or 0.5 for a level that is already filled, empty, or half-filled (containing a single 
nucleon) respectively. In this formulation, the weighting function corresponding to the 
availability of the proton (neutron) level structure of the daughter nucleus, in (3~^ 
depends on the energies of both the daughter nucleon levels e k and the parent nucleon 
levels £j, because their difference is necessary to determine if the decay is energetically 
possible: 

W(e h e k ) = n k Q(Ei -e k ± A NP - m e ) (5) 

where the ± is for /J^-decay, m e is the electron mass, and A^p is the neutron-proton 
mass difference. 

Each composite state is treated as a mixture of quantum states that will depend 
on the model used. Therefore, the mixture in a level i is characterized by a sum of 
eigenstates weighted by the coefficient u>(i;n,l,j) as a function of the quantum numbers 
of n, I, j. For this study, the quantum numbers are for eigenstates of a spherical shell 
model based on a Woods-Saxon potential. The coefficients for a level i are normalized 
to unity: 

^2uj(i;n,l,j) = 1 (6) 

n,l,j 

for each level i. The description and calculation of weights, as well as the determination 
of single-particle states, are further discussed in §2.21 

Transition strength functions are then adjusted to accomodate the single particle 
states in the daughter nucleus. Discrete functions are used to specify the most probable 
location and quantum numbers of an available state. Then, the transition strength for 
a specific level i can be written as a sum of the product of weights for the parent and 
daughter levels, since levels in the daughter nucleus are now assigned. The transition 
strength for favored decays now becomes: 

D n (E,e) = D n (e/,ei) = 5>(i,CM/,0MC,0 (7) 
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where the function A is used to designate the selection rules (for the set of quantum 
numbers ( and £ of the initial state i and final state f respectively) which satisfy the 
transition type Q: 

Af =SUf,3i)S(lf,U)S(T 0J ,T Q)i ^l) 
Act = 5 (j f ,j { ) 5 (l f , k) 5 (r 0J , r 0ji =F 1) 

+ 5 UfiJi ±1 ) s ( l fi l i) 5 (t ,/, T 0:i =f 1) 
Ai,y = 5 (j f Ji ±1)5 (If, k±l)5 (t 0) /, r 0jl =F 1) 
Ai,a = 5 {jf,ji)S(lf,k ± 1) 5 (t ,/, T ,i =F 1) 

+ 5 (if,], ±1)5 (If, U±l)5 (t ,/, r 0il =F 1) 

+ 5 (jf,ji ± 2) 5 (If, k±l)5 (t 0)/ , r 0ji =F 1) (8) 

The strength functions are also normalized to satisfy the necessary sum rules |21j . 
For a single energy level, the Fermi and GT strength functions are normalized to unity: 

/oo 
Df,gt(E, e)dE -> ^ D F ,ar(e f , e<) = 1 (9) 
/ 

while the vector and axial first-forbidden functions are normalized as 121 1: 



J]ZV,a(^) = ^ (10) 

/ 

where R is the nuclear radius. By substituting into Equation [lj the decay rate for f3 T 
decay is now: 

m, 5 c 4 f 

A = g2 ^ 3 ^ 7 \ ^2 n i n k®( £ i - £ k ± ^TVP ~ ^e) 

x EEE^0^;e)Aa(C,e)|Gn| 2 (2Zo + l)/n(^)| (11) 

The variables C, £, and Aq(£,£) represent the quantum numbers associated with the 
states comprising the parent and daughter states respectively and the factor to satisfy 
the necessary selection rules, as mentioned above. The decay order is given by I. The 
summation over Q is simply the sum over all possible transition orders in the decay. Gq 
is the appropriate coupling constant (vector or axial- vector) for the transition. The value 
of C is equal to 1 for allowed transitions and (ft/m e c) 2 for first forbidden transitions, 
the only possibilities considered in this paper. The value E is the transition energy, and 
is simply the difference in parent and daughter energy states plus (minus) the neutron- 
proton mass difference and the electron mass in decay. 

The transition matrix elements for a nucleus are then calculated as: 
\M n (E)\ 2 = \M n (e t -e k )\ 2 

= E i EftEcE € nin*w(i;CMA;;OAi 1 (C,0 1 ) 

The daughter level k is determined by the parent level i, the transition energy E, and 
the transition order Q. The sum over quantum states ( in the i th parent level as well 
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as the sum over states £ in the ] th daughter level are both necessary. Though the 
daughter state is constrained by the transition selection rules for a transition Q, there 
may be several quantum states so allowed, so the sum over states £ is necessary. As 
an example, the discrete transition strength matrix elements for 132 Sn are shown in 
Figures [Q - El One should note the characteristic broad peak corresponding to the 
GT matrix element (as compared to the Fermi matrix element, which has no viable 
transition strengths relevant for (3~ decay), as well as the double peak characteristic of 
the first-order transition elements. 

It can be seen from Equation [11] that excited states are completely accounted for 
in the values n, and n^. This makes the utility of the equation obvious. As well, nuclear 
excitations can take the form of proton or neutron single-particle excitations without 
regard to the type of decay. The inclusion of selection rules is accomplished through 
a specification of quantum numbers of individual states. Furthermore, if all of the 
possible energy states above the Fermi surface are specified, the population of these 
states are assumed to be in thermodynamic equilibrium, and the partition function is 
then determined. Therefore, at temperatures peculiar to the r-process environment (or 
any environment), the values n^ and n^ can be substituted with the probability functions 
for a statistical ensemble in thermodynamic equilibrium. 

2.1.1. Form Factors Besides the calculation of single particle energies, the only 
remaining term in the rate equation is the f-function which, based on the order Al, 
rank Aj, and vector nature of the transition (axial vector or vector) for /^-decay is: 



where Fi,a/v is the Fermi function, which is taken to be the same in all cases; F/ jJ 4/y=F , 
independent of Al. The quantity Sj^a/v is the factor for the Coulomb interaction 
between the outgoing electron wave function and the nucleon wave functions. W is 
defined as the electron kinetic energy E e /(m e c 2 )±l for /3 T -decay, and Wo corresponds 
to the maximum possible electron energy E for a transition. Also, the terms p and q 
are given by p=(W 2 -l) 1 / 2 and q=W -W. The forms of the wave function factors used 
are those of reference [15] . 

2.2. Calculation of Single- Particle States 

2.2.1. Energy Levels of the Single- Particle States The finite range droplet model 
(FRDM) with the Lipkin-Nogami (LN) pairing force is used to calculate single-particle 
levels for the various nuclei. The formula is referred to as a microscopic-macroscopic - 
or "mic-mac" model. The macroscopic part is a smooth function of proton number Z 
and neutron number N based on a finite range liquid drop model. The microcopic part 
is used to calculate individual shell energies and is based on a folded- Yukawa single- 
particle potential. The macroscopic part is described if reference [22]. It is based on 





(13) 
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the overall charge, mass, and the nuclear shape. The parametrized macroscopic part, 
as described in the reference is: 



E mac (Z,N, shape) = 

+AM H Z + AM n N 

+ (-01 + JS 2 - \Ke 2 ) A 



+ 



a 2 + ifP 



9 J 2 A2;B| 
Si 

+a 3 A 1 / s B k 
+a A° 

+ c i 7173-^3 
-c 2 Z 2 A 1 / 3 B r 

2 B W B S 



A 2/3 



-c 5 Z'< 
+/0-4- 



Bi 



-Ca (N 

+W 



+ < 








Z=N=odd 

otherwise 
Z and N odd 
Z odd, N even 
Z even, N odd 
Z and N even 



mass excess of hydrogen and neutron 
volume energy 

surface energy 

curvature energy 
A energy 
Coulomb energy 
volume redistribution energy 
Coulomb exchange correction 
surface redistribution energy 
proton form factor correction 
to Coulomb energy 
charge asymmetry term 

Wigner energy 



average pairing energy 



(14) 



—a P iZ 



2.39 



energy of bound electrons 

The parameters for this expression, as well as their derivations, are described extensively 
in [22], and will be described briefly here. Many of the constants used in this relationship 
are described in Table [3j Other derived quantities are given in the table, as well. The 
standard droplet model defines the shape- dependent quantities [23] as the surface energy 
B s , the surface redistribution energy B w , the curvature energy B k , and the volume 
redistribution energy B r . These are solved analytically by integrating the appropriate 
quantities over the surface of the nuclear volume: 



B s 
B w 
Bk 

B r 



225 
647r 3 rgA 2 

1 

8-KroA 1 / 3 

1575 

647r3r<;A 7 /3 



J s dS 

Is 
Is 



W(r) 



1 

Ri 



1 

R 2 



W(r) 



dS 
dS 
d 3 r 



(15) 



where R\ and R 2 are principal radii of curvature (two principal radii assume 
deformations having no higher than quadrupole terms in the macroscopic model). The 
volume term is specified as W (r) = W(r) — W and: 



W 
W 



I 



V \r~rl\ ( P r ' 

3 IyW( 



d?i 



(16) 
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The quantities B\ and B 3 are the relative generalized surface and Coulomb energies 



respectively for a deformed shape of volume V — ^ (tqA) 3 . These quantities are 
developed to treat nuclei under finite compressibility and diffuseness ad en - The 
generalized forms of these quantities are: 

Bl =^^ff (2-^) e p^d 3 rd 3 r' 

1 87H7-Qa 4 J JV \ a I \r—r'\/a «\ 

i) _ 15 /• /• ■,!'•<' r ( ■ ' ~- r " ~ x 1 - ' -' 1 ' " ' 

• ! 32 7 r 2 r5 J 4 5 /3 J JV \r- 



1 _|_ 1 1^1 ^ p\?~f\/a de 



2 Q*den 



These shape factors, along with their derivatives with respect to the radius are 
solved numerically in the FRDM, and the balance between the Coulomb energy, 
compressibility, and the surface tensions are found based on shape parameters taken 
from the deformation parameters in a folded Yukawa potential. For a spherical geometry, 
Bi and B 3 are equal to 1 with the ranges a and a<j en are equal to - a sharp surface. 
The macroscopic constants in the FRDM have been determined based on a non-linear 
least squares fit requiring about 1000 iterations of the process of choosing parameters, 
recalculating potential surfaces, finding the fit, and re-adjusting parameters. 

The microscopic portion of the model involves the individual shell terms plus the 
effective interaction pairing gaps for protons and neutrons [2H [25] : 

a (18) 

^G p — Z l/3 

in which a constant level density is assumined in the vicinity of the Fermi energy level. 
The value of r mic is determined by a least-squares fit. 

Finally, microscopic corrections to the FRDM are carried out by decoupling the 
shell and pairing potentials as well as the proton and neutron potentials. The single- 
particle potential is then the sum of the central potential V%, which - as stated - is the 
folded Yukawa potential, the spin-obit potential V S , Q _, and the Coulomb central potential 
for protons Vq- Two-particle levels as a function of nucleon number are compared to 
those determined by the Nilson model (with no deformation) for the 132 Sn nucleus in 
Figure H] showing a reasonable agreement between the two. 

For a more complete description of the LN pairing interaction and microscopic 
interaction, the reader is referred to references [22] and to [25] and references therein. 

In the present calculation, single particle levels are determined for both the parent 
and daughter nucleons. That is, transitions between specific levels are calculated for 
a nucleus with j parent nucleons and k daughter nucleons to a nucleus with j-1 parent 
nucleons and k+1 daughter nucleons. Thus, corrections to any shell and pairing energies 
are intrinsically accounted for. 



2.2.2. Quantum Numbers of Single Particle States A knowledge of the quantum 
numbers of each level enables one to calculate the order of a particular j3 transition. 
However, as stated previously, a Fermi Gas model level at a specific energy may 
correspond to several shell model levels in the present formulation, so it may turn 
out to be a mixture of single-particle shell model eigenstates. The basis of states is 
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assumed to be those of a spherical shell model with a Woods-Saxon potential The 
weight for each set of quantum numbers in a level is determined by the proximity of a 
level's energy to that of a basis state. The overlap of the eigenfunctions at a shell level 
determines the strength of their presence in the level (Figure [5]). The method of Nakata 
et al. is used to determine not only the quantum numbers of parent nucleons, but also 
those of daughter nucleons Triangular distribution functions centered on the energy 
of the states are defined: 

' I 1 - ll ^r] < £ " ^ < Q(A) 

Hnjjlc)-- { (19) 
otherwise 



where e n) y is the energy of the eigenstate corresponding to the usual quantum numbers. 
The function is already normalized to N n ^j, which is equal to the degeneracy of a level, 
2j+l. The width q(A) of the basis function of the i th level in a nucleus with \x levels is 
given by: 

q(A) = 6^ (20) 

where d, is the standard level spacing: 

2 

di = — r£° (21) 

1 3 (2i — 2) 1 K J 

where e® is the energy of the Fermi surface of the i th level. Since the values of e are 
allowed to vary smoothly, and scale smoothly with level, Fermi gas levels are appropriate 
to effect this scaling. 

The total sum of the functions in Equation [19] is also defined: 

G(e) = J2Hn,iAz) (22) 

n,l,j 

And G(e) is automatically normalized to the total number of nucleons. To find the 
coefficient u;(i;n,l,J) in Equation [6] for the i th single-particle level, it is necessary to 
normalize the function G(e) such that the overlap in the the region of the level is equal 
to two - the number of nucleons per level: 

Ni J G(e) de = 2 (23) 
where the integration limits are the midpoints between the i th level and adjacent levels: 
s ± = £t + ^ l±1 (24) 

For i=l e~=0. The coefficient for a particular configuration of quantum numbers in the 
i level is then: 

1 f £+ 

^Wn,l,j) = —J H n ^j{e)de (25) 

This method is diagrammed in Figure [5] [9] . The degeneracy of the spherical shell model 
level £ n ,i,j is already accounted for in Equation [TH] with the factor N n) jj. The selection 
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rules in Equation [TT] are then given by the function A^(£;£) and can now be easily 
determined, where Q is the order and rank of the transition (AZ,Aj); ( and £ represent 
the quantum numbers of the initial and final nucleon states respectively, (nj,Zj,jj;n/,Z/J f). 
Thus, the selection rules can be simply stated as: 



where the values of Al and A J depend on the initial and final state quantum numbers. 

Because each state's wave function is a sum of shell-model eigenstates, the smearing 
of single-state quantum numbers over several single-particle states is possible. The net 
result is still an average calculation of /3-decay rates, but now individual particle and 
particle-hole excitations can be simulated by promoting an arbitrary nucleon to any 
available level. Further, the selection rules are now easily calculated using the method 
of the previous section. 

2.2.3. Corrections in Pairing Energy Due to Decay If the highest-level paired neutron 
decays, the other neutron in that level gains an amount of energy roughly equal to the 
pairing energy of that level. As well, paired neutrons in the next lowest level are then 
expected to reconfigure by both losing an amount of energy equal to the pairing energy 
of that level. Thus, the endpoint energy (taken to be positive in this respect) is changed 
by an amount equal to 2b M _i-b^ just from the changes in the parent nucleus. (It is not 
necessary to include the energy of the decaying neutron since this is accounted for in 
the quantity £j — £/.) 

Similarly, in the daughter nucleus, if the highest proton level (in /5~-decay) is 
unpaired, and a neutron decays to that level, then the highest original paired level 
loses its pairing energy for the two protons, while the level to which the neutron decays 
becomes the highest paired level, increasing the total endpoint energy by an amount b M , 
and the total energy available to the ejected electron changes by b At -2b M _i. (Again, the 
pairing energy of the decaying particle is not included for the same reasons mentioned 
above.) 

2.3. Results of (3 -Decay Rate Calculations 

Decay rates and Q- Values have been calculated with the FRDM in previous works [27J. 
The level of accuracy of these calculations is shown in Figures [B] and [7] for about 200 
sample neutron-rich nuclei unstable against j3~ decay for which the decay rates and Q- 
values are experimentally known. In Figure O the value of S is defined as: 



In each figure, most calculated rates fall within an order of magnitude of the true rate. 



A{n h kji] n f , If, j f ) = 5 (l f , k±Al) 5 (j f ,ji ± Aj) 



(26) 




(27) 



2.3.1. Excited State Decays A next logical step is to calculate decay rates for nuclei in 
excited states. The present formulation gives no preferential treatment to nucleons in 
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single-particle excited states. One simply needs to know the values of the excited levels 
and the vacant levels for Equations H] and [5j The levels are already calculated from 
§2.21 Levels can be calculated beyond the maximum filled level in the ground state, 
although the simplicity of the model can render these values slightly different from their 
experimental counterparts. 

In stellar environments, the probability of finding a nucleus in an excited state 
at a given temperature can be obtained from the partition function by knowing the 
spin and energy of the excited state. Since each particle level holds only two nucleons, 
the degeneracy of each level is two, and the quantum numbers are given by the level's 
proximity to those of the spherical shell model (as described in § §2.2. 2D . The difference 
in energy: 

m fJ, 

AE = E <A - E <A (28) 

*JV,Z=1 »JV,Z=1 

where values with an * are those for an excited-state nucleus, while those with a are 
those corresponding to a ground-state nucleus. The summation is taken over neutron 
and proton levels. The value m is the highest bound state level. The first term represents 
the total energy of the excited nucleus, and the second term is the total energy of the 
ground-state nucleus. The average decay rate of an isotope in a stellar environment is 
then the weighted sum over energy states [28J: 

A = ^P(A£)A(A£) (29) 

AE 

where P(AE) is calculated using the partition function. 

Excited states and their decay rates have been calculated for several nuclei. A 
sample of the decay rates as a function of these states are shown for four nuclei in 
Figure [HJ This figure shows the dependence of decay rates on excitation for an even-even 
nucleus ( 150 Ba), an odd-N nucleus ( 149 Ba), an odd-Z nucleus ( 147 Cs), and an odd-odd 
nucleus ( 150 La). One first notes the multiplicity of levels resulting from the splitting of 
states with spins higher than 1/2, a direct result of using levels with a degeneracy of 
two. 

One also notes the apparent clustering of decay rates into bands. In the case of 
150 Ba, one notes three major bands at ^0.75, 1, and 2 where Ao is the ground- 
state decay rate. This is due to the fact that single-particle levels are dominated by 
specific shell-model eigenstates in certain regions, as expected. In the case of 150 Ba, the 
first several single-particle neutron excitations above the ground-state contain large 
admixtures of the 213/2 state, while those near the ground-state Fermi surface are 
predominantly in the hg/2 state. In the case of transitions in this region, the decay 
order of the particle-hole configurations goes from GT to first-forbidden, thus reducing 
the decay rate. Single-particle excitations in certain energy regions are expected to 
have similar decay strengths. The band at ~2 is due to proton ng 9 /2vrg 7 72 particle- 
hole configurations which do not change the transition order, but change the decay 
Q- value by opening up a hole at lower energies. The band at A/A ~1 corresponds to 
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^13/2^11/2 configurations. Since this has very little impact on the highest-energy GT 
(or first-forbidden) transitions (to the ng-j/2 or 7rg 9 /2 state), the overall transition rate 
is changed very little. It is expected that higher excitations will result in other band 
structures. However, since astrophysically interesting temperatures are below 1 MeV, 
the contribution to average decay rates from states much above this is small, so only 
lower excitations were studied. 

Similarly, the band structures in the other nuclei shown can be explained in a similar 
manner. The band structures of other nuclei vary based on their individual structures. 
For example, 150 La has a band at ^ < 1, but no band at ^ > 1 (for the energies of 
interest). This nucleus has a higher deformation than 150 Ba, as well as an unpaired 
neutron and an unpaired proton. The lowest excitation energies are accomplished via 
neutron 1^/19/2^13/2 particle-hole configurations, but single-proton excitations are due 
primarily to the unpaired proton, which has little bearing on the available hole, but 
does account slightly for the vertical band structure in the figure; the rate increases 
with the transition energy for a single type of transition. 

Of course, excitations can change the strength functions, as the order of a particular 
transition can be altered. In the case of 150 Ba, the GT strength function is shown for the 
ground state, a transition to the region dominated by the Ttg^/2^9i/2 state, and an excited 
state dominated by the ^9/2^13/2 state in Figured! As expected, the neutron excitation 
lowers the overall GT transition strength, particularly at the highest transition energies, 
where the transition rate is highest. Conversely, the proton excitation increases the 
overall rate slightly at the lowest energies, while also increasing the Q-value. It should 
be noted in this case that, although the rate is doubled and the GT strength function for 
the proton excitation looks very similar to that of the ground-state, the average lifetime 
is still relatively low (~300ms), and the result is a small absolute increase in lifetime (to 
about 315ms). Overall, there are similar small changes in the first-forbidden strength 
functions, shown in Figure [TUl 

Given these figures, one cannot definitively state that excited-state nuclei decay at 
a higher rate than ground-state nuclei, as the transition order shift becomes important. 
The average decay rates of the nuclei in Figure [8] are plotted as a function of temperature 
in Figure HH calculated using Equation [29J The decay rate of 132 Sn as a function of 
temperature is also plotted. While each nucleus behaves differently, it seems clear that 
the doubly-magic 132 Sn does exhibit a sharp increase in decay rate with temperature 
at higher temperature, most likely due to the fact that 132 Sn has closed proton and 
neutron shells. All single-particle excitations in 132 Sn are quite high in energy and open 
transitions that were previously Pauli-blocked. The transition strength functions for the 
first calculated excited state of 132 Sn are shown along with the ground-state functions 
in Figures [Ul3] indicating greatly increased GT strengths even at high transition energy. 

Finally, the effect of temperature on /3-decay rates is evaluated as a function of 
the neutron-richness of the nucleus in Figure [131 In this figure, relative decay rates are 
plotted as a function of temperature for various Z on the same isobar A=162. For this 
isobaric chain, there may be a slight dependence on the stability of the nucleus. While 
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more study is warranted, this may be of interest in suggesting that GT transitions of 
the highest level neutrons in the very neutron-rich nuclei along the r-process path are 
most likely not Pauli blocked, as they might be for the nuclei closer to stability. In 
this case, single-proton excitations will occupy states that a neutron would otherwise 
decay to. Nuclei closer to stability, on the other hand, have truncated GT strength 
functions. Thus, proton or neutron excitations are more likely to open transitions that 
were originally Pauli-blocked. 

3. Effects of Excited-State Decays on the r-Process 

Relative changes in decay rates along the r-process path are shown in Figure [12] for 
the line of constant neutron separation energy S n ~ 2.5 MeV for three different 
temperatures. The effect is to increase slightly the relative ratios of the ground-state 
rates at lower mass. However, this effect is only slight (<10%) in temperature regimes 
relevant to the r-process. A network calculation is necessary to evaluate fully the 
magnitude of the effect on the r-process. The increase in the 132 Sn rate may be enough 
to make a significant difference in the final abundances, but this is not clear yet. 

3.1. The Explosion Model 

The r-process environment used in the present study is constructed in the semi-analytic 
model of the nutrino-driven winds of reference [31] . The mass flow out of the supernova 
hot-bubble region is treated as a spherically symmetric, steady flow around a neutron 
star of mass M and radius R in the Schwarzchild geometry, which is described [30] by 
the following sets of equations: 



energy, u is the radial component of the four-velocity, and q is the net heating rate; q is 
positive just above the surface of the neutron star, decreasing roughly exponentially with 
radius. The total energy density pr is related to the baryon density by pr = Pb (1 + £)• 
In the above equations, the constants h, c, ks, and G are equal to unity. In a stationary 
rest frame, the radial velocity of the expanding bubble is the product of the Lorentz 
factor and the radial component of the four-velocity: 



The net heating rate is determined by the sum of neutrino heating and cooling from 
the following reactions and their inverses: 

v e + n <-> e~ + p 




(30) 



where M is the mass outflow rate, p& is the baryon density, e is the specific internal 




(31) 



Excited- State (3 -Decays and the r-Process 



15 



V e + P e + 72 

Ui + e~ <-> Vi + e" 
Vi + e + <-> z/j + e + 
z/j + e~ <-> z/j + e" 
z/j + e + «-> z/j + e + 

z/j + z/j <-> e~ + e + (32) 

where the index z corresponds to each neutrino flavor. The Newtonian forms of these 
heating and cooling rates are described in detail in references j6] and [29]. Since the 
mass flow is influenced by the properties of the proto-neutron star, the heating rate is 
strongly dependent on the neutrino luminosities and energies. Transformation of the 
solid angle from the reference frame of the neutron star to the r-process site leads to the 
correct form of the heating rate in the Schwarzchild geometry |31j . In this calculation, 
the neutrino luminosity L Ui is assumed to be the same for all neutrino flavors and the 
RMS energies are used; the energies for e Ue , e Pe , s u ^ Tl and e„ are 12, 22, 34, and 34 
MeV respectively. 

Equations [30] are supplemented with the equations of state in the high-temperature 
limit taken from Qian and Woosley (1996): 

p = — T 4 + -^_T 

180 m N 

e= ^ + lZ. (33) 
60 2m N K ' 

where is the nucleon rest mass. 

Using Equations [30] and [33] the evolution of the hot bubble is followed. The 
adjustable parameters for solving the equations are M, R, and M. The neutrino 
luminosity is taken to be 5 - 7xl0 51 erg s -1 . The equations are solved implicitly and the 
result is the velocity and thermodynamic quantities as a function of r (and hence, t). Any 
trajectory (velocity, temperature, and density as a function of radius) can be specified 
with the parameter set, and the entropy and dynamic timescale - parameters important 
to the r-process - are direct results of the evolution of the r-process environment. In this 
particular calculation, the approximation of a steady flow is used |31j . which allows an 
analysis of the rarified region about the supernova core without requiring a tremendous 
amount of computational power. 

For this model, a typical outer boundary temperature of T=0.1 MeV behind 
the shock at r^l0,000 km, consistent with the theoretical result of the benchmark 
hydrodynamic simulations of the delayed-explosion of Type II supernovae [6] is used. 
Subsequent extended studies of the neutrino- driven wind [32] [33] have used the outer 
boundary condition that the mass ejection rate M (Equation [301) is taken to be 99% of 
its critical value (that is, the value for which wind velocity is supersonic) . We therefore 
adopt the same outer boundary condition as that of ref. [32] with the core radius 
held constant at 10 km in the present calculations is used. The only two remaining 
hydrodynamical parameters are the core mass and the neutrino luminosity. From this, 
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the dynamic timescale r dyn} temperature, and density at any radius can be determined. 
From the temperature and density, the entropy is determined. 

The nuclear reaction network used in this calculation [3U [35] is coupled to the 
output of the aforementioned hydrodynamic code. As shown in Figure HH the network 
consists of about 3700 nuclei with Z<93 and includes all nuclei between the most proton- 
rich stable nuclei and the neutron drip line [31]. Possible reactions in this network are 
(n,7), (p/y), (a, 7), (p,n), (a,p), (a,n), /3-decay, /3-delayed neutron emission, electon 
capture, and neutrino neutral- and charged- current interactions and their inverses 
[31] ■ These reactions are summarized in Figure [T5j The effects of excited-state 
nuclei are included via the scaling of decay rates, as discussed below. Nucleosynthesis 
calculations begin at T 9 =9, a sufficiently high temperature that a state of nuclear 
statistical equilibrium (NSE) exists, defined by the balance of strong and electromagnetic 
interactions. At this temperature matter is mostly in the form of free nucleons. The 
initial electron fraction Y e is a parameter of the reaction network. 

In order to maintain the consistency in the mass formula, the same formula used 
to calculate shell energies for /3-decay rate calculations (i.e., the single-particle formula 
of TUYY p6j)was used to calculate neutron separation energies, a-particle separation 
energies, /3-delayed neutron emission probabilities and ground-state /3-decay rates for the 
nuclei shown in Figure [HI /3-decay rates as a function of temperature were determined 
for nuclei with neutron separation energies less than about 2.75 MeV. (It is not very 
useful to calculate decay rates as a function of temperature for nuclei with S n >2.75 MeV 
because, the r-process path does not pass through these nuclei until well after freezeout.) 
Of course, the neutron separation energies and capture rates are also affected by nuclear 
single-particle excitations. This becomes important especially along the r-process path 
where neutron separation energies are low. One may expect an increase in both the 
neutron emission rate as well as the capture rate for excited-state nuclei though capture 
and photoemission are in equilibrium in the classical r-process. However, single-neutron 
excitations cannot go above the neutron separation energy. For this reason, single- 
neutron excitations are kept realistically low (though multiple-particle excitations may 
be higher). Indeed, even at r-process temperatures, one does not expect a significant 
population of states above the neutron separation energy. 

While the network calculation proceeds to Z=100, the mass formula used in this 
work has shell energies up to N=157. Therefore, neutron-rich nuclei with N>157 cannot 
be quantified, which effectively limits proton number Z to about 93. However, this is 
well past the A=195 peak in the r-process abundance distribution, and the effects of 
not including excited state rates for these very heavy nuclei have little effect on the 
abundances of the lighter nuclei. Also, it was found that the abundances of fissile nuclei 
along the path are very low at freezeout, so that the inclusion of fission cycling was 
ignored in this r-process model. 
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3.2. Excited State (3 -Decays in the Reaction Network 

The inclusion of excited state /3-decays into the reaction network can be computationally 
simple. Functional forms of decay rates with temperature were used. This allows a 
fairly accurate evaluation of the calculated rates of individual nuclei for temperatures 
between and 6 billion K, while still maintaining the speed of the network calculation. 
Sixteen regions of the isotopic chart were found for which decay rates of nuclei in these 
regions have similar dependences on temperature. These regions are shown in Table HI 
Nuclei represented in these figures have neutron separation energies S n ranging between 
and about 2.75 MeV. This sort of parametrization allows for a very reasonable first 
approximation of the change in decay rates over a large region of the isotopic chart. 

The dependence on temperature for the above mentioned regions can be conveyed 
with a set of parameters fit to the formula: 

+ aexp(/3T 9 + 7 T 9 2 ) (34) 

where A is the /3-decay rate at temperature Tg (defined to be the temperature in 10 9 
Kelvin), and the ground-state decay rate is A . 

The ratio in Equation [M] must be unity at Tg=0, so several of the parameters are 
dependent: 

B= {A -I) (exp (-CD) + 1) (35) 

so the parameter set is reduced to eight in number. Further, in using Equation [34] to fit 
the dependences of rates to temperature, two terms were necessary in the temperature 
region of interest. In actuality, only one of the two terms is present at a time. So if E 
is non-zero, then a is zero and vice- versa. Also, a is limited to one of two values - 
or 1 - leaving seven free parameters. The parameters for each of the regions mentioned 
above, along with the nuclei included in each region, are given in Table |U 

Two of the parameters in Table H] are functions of the nuclear proton number Z. 
These are the values of A for region 8 and the value of 7 for region 7a, listed in the 
table as A(Z) and j(Z) respectively. This is a convenient way of representing these 
parameters over a very large isotopic region, as the functional form of the temperature 
dependence of decay rate changes slowly with Z. These parameters are given by: 

A(Z) = 100.31 - 1.38Z 

%Z) = 0.185 - 0.003Z (36) 

No odd-even effects are observed because Equation [34] is a ratio of rates, so any local 
odd-even effects would be minimized in the division. Variations in the parameters as a 
function of mass seem to be explainable in terms of local shell effects. 

3.2.1. Uncertainty Associated With the Use of the Spherical Shell Model It is worth 
mentioning the use and applicability of the spherical shell model in this formulation. 



^ = E(FT 9 + G) 



A 



B 



exp (C(T 9 -£>)) + ! 
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One should not only consider effects on decay rates, but on the overall r-process. While 
it is difficult to compare to a more elaborate model in this work without utilizing such a 
model, a qualitative assesment is discussed here. Two things which are not necessarily 
mutually exclusive must be considered. The first is the accuracy of the decay rates 
used. The second is the applicability of the spherical shell model to the nuclei studied. 
In particular, one can estimate effects on individual transitions when one considers 
a Nilsson shell model with deformation. For this reason, nuclei which are presumed 
undeformed are compared to those which are considered to lie in regions of deformity 
in the isotopic chart (i.e., between closed shells) are discussed. 

In estimating the accuracy of the calculated decay rates, one considers uncertainties 
in calculations for ground-state nuclei; the uncertainty in the decay rates may be 
estimated using Figure EJ The rms error in S (Equation 1271) for all nuclei studied 
in this work with known f}~ decay rates has been calculated. The rms value S rms =(S 2 )^ 
has for all nuclei, even-even nuclei, odd-N nuclei, odd-Z nuclei, and odd-odd nuclei to be 
0.767, 0.569, 0.805, 0.745, and 0.872 respectively for the ground state nuclei, indicating 
reasonable agreement as shown in Figure [6] for at least the nuclei with known decay 
rates including those in regions of presumed deformity. The accuracy improves with 
rate as the presumed level density increases. For nuclei with A >0.1 s _1 the value of 
S rms is 0.579 and decreases further to 0.467 for nuclei with A > 1 s _1 , indicating the 
improvement in the accuracy of this formulation for the r-process, in which the decay 
rates are larger. The agreement is also equally good for the few known neutron-rich 
isomers, though some caution is warranted here as the number of known isomeric states 
is small. As many of the nuclei studied were extremely neutron rich, the uncertainty 
in calculations is expected to be reasonable for those nuclei along the r-process path. 
In fact, the uncertainty is expected to decrease with single-particle level density. While 
reasonable, these errors may likely be reduced with a more realistic model. 

In the case of non-deformed nuclei, the spherical shell model is believed to be a 
reasonable estimate of the decay rates. This is extremely important in the case of the 
r-process as the waiting points are localized about areas of lowest deformation; most 
of the r-process abundance is confined to these regions of the isotopic chart. Large 
uncertainties in the decay rates through the waiting points can result in drastic changes 
in the final distribution. For example, consider the A~130 abundance peak. If the 
decays rates of the r-process nuclei associated with this peak were much higher, then an 
increased flow through this peak would result, resulting in a much lower abundance of 
nuclei in this peak and a very large enhancement of the rare earth nuclei (130<A<190). 

On the other hand, the portion of the r-process path thought to be responsible for 
the production of the rare-earth elements passes through a region of the isotopic chart 
associated with possibly large deformations. The overall effects on the final r-process 
abundance distribution due to changes in the decay rates of these nuclei is expected to be 
small as their contribution to the total abundance is small. In the case of these possibly 
deformed nuclei, single-particle levels in the Nilsson model may shift by as much as 0.5 
MeV from those of the spherical shell model in cases of extreme deformation. However, 
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single-particle levels are shifted in both the positive direction and the negative direction 
for protons and neutrons. Qualitatively, the net effect is to break degeneracies in closed 
shells while "spreading out" high localized single-particle level densities associated with 
closed shells. 

Consider, for example, Figures [9] and [10] which show the calculated GT and first- 
forbidden strength functions for the 150 Ba nucleus for various excitations. If a more 
realistic structure model was used, the net result would be a spreading of the widths of 
the peaks of these functions as both the daughter and parent states would undergo shifts 
in both the positive and negative energies. While some of the strength function would 
shift to higher transition energies (more negative energies on the figure) resulting in 
contributions to higher decay rates, some would also shift to positive energies resulting 
in contributions to lower decay rates (and possibly - though unlikely - to states that 
cannot decay). By shifting the transition energy by an extreme amount of 0.5 MeV for 
a GT transition in Equation [TTJ it is estimated that the transition rates may change by 
as much as a factor of two. However, it must be carefully noted that this is for a shift 
of the entire GT transition strength in one direction for an extreme amount for only 
the most deformed nuclei. Thus, this factor is an upper limit for very deformed nuclei, 
and likely the error due to using a more accurate model is found using the rms values 
in S as discussed above. Certainly future work may concentrate on more realistic model 
calculations, though it will be seen from the next section that this may not be necessary 
as the results of the r-process calculations are more heavily dependent on uncertainties 
in the hydrodynamic conditions. 

4. Results of the Network Calculation 

The results from several hydrodynamic parameter sets, as well as electron fraction 
parameter values Y e , were examined. These parameter sets are shown in Table [5] For 
each parameter set, the core mass in solar masses, core radius, neutrino luminosity, 
initial electron fraction, and whether or not /3-delayed neutron emission is included 
are listed. Using these parameters and the calculations of reference [31] the dynamic 
timescale and the entropy in the expansion are constrained. Though still in agreement 
with current predictions, the dynamic timescales in these calculations are shorter than 
average. However, the entropy is lower, and no artificial increase in the entropy (as is 
often assumed) was required [6]. 

Each simulation is run until several seconds beyond freezeout. While this time is 
sufficient to gauge the gross features of the r-process abundance distribution, a longer 
simulation may have resulted in more post-processing, allowing for smoother abundance 
distributions. Figure [16] shows that model A underproduces the A~195 peak by a large 
amount. This is due primarily to the fact that the lower entropy in model A results in a 
very low neutron-to-seed ratio. Models B and C were chosen as intermediate points in 
the entropy-timescale phase space. Both produce a more acceptable r-process abundance 
distribution, although the A~195 peak is still underproduced. For comparison, the solar 
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r-process abundance distribution is displayed in the figure. One notes some residual 
even-odd effects in the calculated distribution as the effects of smoothing may not be 
complete, though the gross features of the distribution are noted. 

Results from models B and C are displayed in Figure [T71 for both the hot (i.e., 
including excited-state decays, solid line) and cold (i.e., not including excited state 
decays, dashed line) models. The A~195 peak is most profoundly affected, along with 
the nuclei just below this mass region. From the relationship between the decay rates 
and temperature (Equation [3H and Table Hj), it can be seen that decay rates of the 
nuclei in the region just below the A~195 peak (and - to a lesser extent - the region 
just below the A~130 peak) are quite sensitive to changes in temperature even at low 
temperatures. This is expected due to the high level densities of these nuclei (lying just 
below the N=126 and N=82 closed shells), as discussed in £j2l Shell quenching has not 
been included in this calculation, though the effect is noted, and no conclusions can be 
drawn from this study to evaluate the effects of quenced closed shells far from stability. 
Other rates, however, are not as sensitive to temperature changes at low temperatures 
and, as the r-process progresses, these rates would drop to their ground-state values 
before those of the nuclei in the regions below the abundance peaks. The decay rates 
of nuclei in this region would increase relative to those of nuclei in other regions of 
the path, selectively depleting the abundances of nuclei in this region. This effect is 
displayed in Figure Q2I in which nuclei in these two regions have large changes in decay 
rate as Tg increases a small amount. 

The lowering of the electron fraction in models F and G is a physically acceptable 
assumption given that the electron neutrinos and anti-neutrinos captured on nucleons 
before the beginning of the a-process can alter the neutron excess rj via neutrino 
interactions on the free nucleons in NSE. In a steady-state, the electron fraction is 



Given the neutrino energies assumed in the nucleosynthesis code, then a value as low as 
Y" e =0.35 may not be unexpected, though higher values have also been used. 

The abundance distribution results of models F and G are shown in Figure [181 
The A~195 peaks are more pronounced - closer to acceptable values - in both cases as 
compared to models B and C. Even more, abundances of nuclei with heavier masses are 
also increased, more closely matching the solar system abundance distribution in this 
region. As in previous models, the effect of excited state /3-decays is still quite evident in 
both models, especially for the heavy nuclei with A>195. The ratios of peak abundances 
for these models are shown in Table El Both models have peak abundance ratios close 
to that of the solar system. Also, in both models, the A~195 peak is shifted closer 
in position to the solar system abundance peak. The slight decrease in Y e , while still 
maintaining physically realistic values, combined with the use of excited-state /9-decays, 
produces an r-process calculation with results that match the solar system distribution, 



[37j 




(37) 
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with the exception of nuclei in the A~180 region, which seem to exhibit a deficiency in 
abundance. 

Figures H7] and [[8] suggest that the inclusion of excited-state /3-decays results in a 
net shift of abundance from the region of nuclei just below the A~195 abundance peak 
into this abundance peak. Certainly, one notes that the more rapid decay rates in the 
region just below the abundance peaks results in a faster progression of the r-process 
through these nuclei. However, one should also note that the integrated abundances for 
nuclei with 180<A<200 is higher in the hot model by than in the cold model by 29%, 
indicating that overall abundance is increased in both the abundance peak and the region 
below this peak. This suggests that the increase in abundance in the A~195 abundance 
peak is only partially due to an increase in the decay rates of the nuclei associated 
with the region just below this peak. One must note that the overall abundance in 
the 180<A<200 region increases due to a flow of abundance into this region from the 
lower masses. This increased abundance can come from minute changes in abundance 
in the A~130 peak and from the mass region just below this peak. In the semi-log 
representation of Figures [T7J and [181 it is seen that even small fractional changes in the 
abundance of the A~130 peak - with nearly an order of magnitude more abundance 
than the A~195 peak - can result in sizable changes in the 180<A<200 mass region. A 
similar comparison can be made between the abundance region just below the A~130 
peak and the A~195 peak. 

This brings up an interesting point regarding the net effect of rate increases with 
temperature in the r-process. Not surprisingly, figure [12] shows that decay rates for the 
closed-shell nuclei are only affected at high temperatures (corresponding to early stages 
in the r-process). Very early in the r-process, flow through the A~130 mass region 
increases relative to the remaining nulcei along the r-process path, populating the rare 
earth region. However, as temperature drops before freezeout, the relative flow through 
this mass region decreases, and the relative flow through the rare earth region is still 
higher, populating nuclei in the A~195 mass region. As the A~195 mass region is not 
populated until later (and cooler) in the r-process, the effects of excited-state decays are 
minimal. The net result is that as the r-process progresses to higher mass, the relative 
decay rates drop on average first for the A~130 abundance peak and then for the rare 
earth region. Abundances of the r-process progenitor nuclei are affected early on by 
increased rates of the A~130 nuclei and later on by those of the rare earth nuclei. Of 
course, the net effect is dependent on the passage of the r-process path through many 
nuclei, so the network calculation is employed as a useful tool. 

Models D and E are identical to models C and B respectively, except for the fact 
that /3-delayed neutron emission reactions are included in models D and E. The final 
freezeout abundance for model D is compared to that of model C, and the final freezeout 
abundance of model E is compared to that of B in Figure [19] The distribution seems to 
be shifted slightly to heavier mass when /3-delayed neutron emission is included. While 
this may be surprising, just after freezeout, as nuclei begin to /3-decay back to stability, 
/3-delayed neutrons become available for capture by all nuclei. The two-neutron emission 
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probability as shown in Figure |2"U1 is higher in the 125<A<180 region than in the region 
surrounding the A~195 peak, while the single-neutron emission probability is slightly 
higher in the A~195 region, resulting in a larger net loss of neutrons in the lower mass 
nuclei. Furthermore, the heavier nuclei tend to have higher neutron capture cross- 
sections [3BJ- The overall result is that the masses of the heavy nuclei are increased 
by a few mass units. Figure [19] shows that the mass shift of the abundance peak is 
roughly two mass units, indicating that nuclei in the A~195 mass region have captured 
at least two net neutrons while decaying back to stability. This necessitates further 
study regarding the dynamical nature of the r- process after freezeout [39] ; the r-process 
path continues to evolve to its final distribution even after the neutron abundance has 
dropped several orders of magnitude. 

5. Conclusion 

This work provides a study of the effects of excited state /3-decays on the r-process. A 
preliminary method was used to evaluate the possible effects of /?-decay rates of excited- 
state nuclei. Though the accuracy of the model is limited by the knowledge of single 
particle levels, an approximate treatment allows one to gauge the magnitude of effects 
on the r-process and provide impetus for further study. An empirical calculation was 
employed to find single-particle levels, and quantum numbers were deduced based on the 
level proximity to those of the spherical shell model. Although minor effects were found 
to result from inclusion of the excited state decays, there are also other possible effects 
that the inclusion of excited state nuclei may have on the r-process. One can imagine 
that if the excitation is due to the promotion of neutrons to higher-lying single-particle 
orbitals, then the photoneutron Q-value will decrease, and the (7,11) reaction rate may 
increase. Thus the effect of an increased rate might be to shift the r-process path closer 
to /3-stability. 

All models used in this calculation do a reasonable qualitative job of reproducing 
the solar abundance distribution for the mass region 80 < A < 130. The fact that the 
abundance at low mass is roughly independent of the type of model used (hot or cold) 
is an indicator that the nuclei in this region have decay rates not as heavily dependent 
on temperature as some of the higher mass nuclei. In the more massive nuclei, the 
sensitivity of the decay rates on temperature resulted in a more pronounced shift in the 
path as the r-process evolves through freezeout. 

As mentioned previously, the dynamical treatment of the r-process is an important 
factor here in that the path continues to evolve even during freezeout, a result of /3- 
delayed neutron emission. With the mass formula used in this evaluation, it was found 
that the low mass nuclei have a higher probability of emitting two neutrons during f3- 
decay than the higher mass nuclei, which have a higher probability of emitting a single 
neutron during /3-decay. These available neutrons are recaptured, with the cross section 
roughly increasing with mass. The net result is a slight shift in the A~195 abundance 
peak. 
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Despite the ability to predict a more reasonable abundance of the A~190 nuclei, 
there is still a discrepancy between the predicted abundance distribution in the rare 
earth region and that of the solar system. The abundances of the rare earth elements 
(A~165) are predicted to be lower in abundance by about an order of magnitude than 
the A~130 peak, in fair agreement with that of the solar system. This corresponds to 
the argument of the authors of reference [39J, who state that the rare-earth region is a 
robust feature of any dynamical calculation including post-production of the r-process 
progenitors. However, nuclei with 130<A<160 are overproduced slightly, removing the 
effect of the rare earth region being manifest as a peak, hence the appearance of the 
abundance distributions in the figure. 

It is obvious that no shell quenching has been employed in this preliminary model, 
as can be seen by the dip in the A~180 abundance. While the abundance of the A~180 
nuclei relative to the A~130 peak is similar to that of the solar system, the width 
of this dip is greater than that of the solar system abundance distribution. It has 
been mentioned [10] that the underproduction of nuclei in this region might vanish if 
the quenching of shell closures for the very neutron-rich nuclei in this mass region is 
properly included. 

Because of the large number of nuclei involved, the improved semi-gross theory 
has been used to globally calculate decay rates. It is understood that the initial model 
represents a first attempt to gauge the effects of nuclear /9-decays in the hot environment 
of the r-process, and further study is warranted. In particular, a more accurate global 
calculation of decay rates is desired. Currently, the measurements of /3-decay rates are 
limited to either ground-state nuclei or long-lived isomers. However, with the advent 
of large neutron flux devices [41], the measurement of the GT strength functions of 
nuclei in excited states from may become feasible within the next decade, allowing for 
experimental confirmation of the transition strengths of excited-state nuclei. 
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Table 1. Terms Used in Gross Theory Presented in This Paper 


Term 


Description 


Equations 


AT 

Mi 


Number of Parent Nucleon States 


1 


w 


Weighting Function Describing 


1,5 




Availability of Daughter Nucleon States 




£ 


Energy of a Nucleon State 


1-4,4,17,18,20 


Do 


Transition Probability Function 


1,7,9,10 


£a 


Lowest Energy of a Parent Nucleon 


1 




Energetically Able to Decay 




£l 


Highest Energy of an occupied Parent Nucleon State 


1-3 


n 3 


Number of Daughter Nucleons Collected at the 


3 




Highest Daughter Level Due to the Pairing Force 




2A 


Width of the Pairing Gap 


2,3 



Table 2. Terms Used in the Single-Particle Rate Calculations Presented in this 



Term 


Description 


Equations 


n. 


Number of Parent or Daughter Nucleons in Level i 


4,5,11,12,21 


CO 


Spherical Shell Model Eigenstate Mixing 
Coefficient for Mapping to Leves 
in the FRDM 


6,7,11,12,20 


An 


Function for Describing Selection 
Rules of a Transition f2 


7,8,11,12,21 


gfi 


Multiplicity Spin- Averaging Factor for 
Transition Type 


11 




£ k 


Discrete Uncorrected Fermi Energy of 
the k th Level 


15,16 


dk 


Standard Level Density About the k th Level 


15,16 


H(e) 


Triangular Distribution Function Indicating 

the Strength of the Spherical Shell Model Eigenstates 


14,17,20 


q(A) 


Width of the Triangular Function H(e) 


14,15 


G(s) 


Sum of H(e) Over All Eigenstates 


17,18 


e ± 


Integration Limits About G(e) and 
H(e) for Normalizing and Determining uj. 


18-20 
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Table 3. Parameters Used in the FRDM Nuclear Mass Model 



Parameter Description Value Units 



AM H 


Hydrogen atom mass excess 


7.289 


McV 


AM„ 


Neutron mass excess 


8.071 


MeV 


O-el 


Electronic binding constant 


1.433 x 1CT 5 


MeV 


K 


Nuclear compressibility constant 


240 


MeV 


W 


Wigner constant 


30 


MeV 


a<3 


^4° constant 


0.0 


MeV 


ax 


Volume-energy constant 


16.247 


MeV 


ai 


Surface energy constant 


22.92 


MeV 


a-i 


Curvature energy constant 





MeV 


J 


Symmetry energy constant 


32.73 


MeV 





Effective surface-stiffness constant 


29.21 


MeV 




Charge asymmetry constant 


0.436 


MeV 


I 


Relative neutron excess 


N-Z 


_ 


A„ 


Average neutron pairing gap 


r ma cB s 
N 1 / 3 


MeV 


A, 


Average proton pairing gap 


z 1 / 3 


MeV 




Average neutron-proton interaction energy 


h 

B A 2 / 3 


MeV 


Cl 


Coulomb energy constant 


3.k. 

5 r 


MeV 


C2 


Volume redistribution energy constant 


336 V .7 T K ) L l 


MeV 


C 4 


Coulomb exchange constant 


l(^) 2/3 -l 


McV 


C5 


Surface redistribution energy constant 


1 r 2 


MeV 


/o 


Form factor correction constant 


1 /145\ r p e2 
8 V 48 ) 


MeV 


ro 


Nuclear radius constant 


1.16 


fm 


r p 


Proton rms radius 


0.80 


fm 


Tmac 


Average pairing gap constant 


4.80 


MeV 


h 


Neutron-proton interaction constant 


6.6 


MeV 


Q<den 


Diffusencss of Yukawa function 


0.70 


fm 




LN pairing constant 


3.2 


MeV 
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Table 4. Parameters Used for Fitting Decay Temperature Dependence 

Parameter 

Region Label 







A 






E 


p 


(y. 


ft 

p 


7 


67<N<75 


1 













i 





0.03 


76<N<77 


2 


2 


2 


2.5 


0.3 


1 









N=78 


3 


3 


9 


0.25 


0.095 


1 





- 


- 


N=79 


4 


1.25 


5 


0.5 


0.2 


1 









80<N<81 


5 













1 





0.018 


N=82 


6 


_ 


_ 


_ 





_ 


1 


0.032 


0.045 


39<Z<42 


7a 


- 


- 


- 





- 


1 





1{Z) 


47<Z<51 




















43<Z<46 


7b 


9 


1 


2.5 


0.2 


1 


o 






82<N<126 52<Z<54 


7c 


12 


1 


4 


0.28 


1 









Z=55 


7d 













1 


0.08 


0.025 


Z>56 


7e 


12 


8 


0.5 


0.04 


1 









N=126 


8 


A 


4 


1.2 


0.3 


1 









60<Z<63 


9a 


5 


3 


0.75 


0.25 


1 









64<Z<66 


9b 


20 


0.88 


7.3 


1 












N>126 68<Z<71 




















Z=67 


9c 


1.5 


5 


0.25 


0.4 


1 









72<Z<74 


9d 













1 


0.3 


0.025 



Table 5. Parameters Used in the Supernova Model 



Model 


M/M© 


R(km) 


L„ (10 52 erg s- 1 ) 


Y e 


CM? 


T~dyn (s) 


S (k) 


A 


1.7 


10 


0.7 


0.4 


N 


0.017 


110 


B 


2 


10 


0.7 


0.4 


N 


0.014 


151 


C 


2 


10 


0.5 


0.4 


N 


0.021 


162 


D 


2 


10 


0.5 


0.4 


Y 


0.021 


162 


E 


2 


10 


0.7 


0.4 


Y 


0.014 


151 


F 


2 


10 


0.5 


0.37 


Y 


0.021 


162 


G 


2 


10 


0.5 


0.35 


Y 


0.021 


162 



Table 6. Relative Abundances of the A~195 and A~130 Peaks in the r-Process 
Calculations (Hot Models) 



Model 


A m aa:,195 


A-max, 130 


Y(A ma:! ; i i95)/Y(A ma:! ; i i3o) 


^A-195/^yl-130 


B 


198 


128 


0.04 


0.05 


C 


196 


127 


0.07 


0.09 


D 


198 


128 


0.13 


0.12 


E 


198 


129 


0.12 


0.13 


F 


194 


127 


0.22 


0.23 


G 


194 


127 


0.27 


0.27 


Solar 


195 


130 


0.27 


0.29 




i 



_ m 



> 



i 



in 



oo 
O 



O 



^1- 

o 



CM 

o 



d 



^1- 

d 



d 



in 



IT) 

d 



Figure 1. Transition strength function \Mct(E)\ 2 assuming discrete levels for 132 Sn. 
The top graph is the discrete transition strength function. The middle graph is a 
histogram of the discrete function with 0.5 MeV bins, and the bottom graph is the 
integrated function. In the lower two graphs, the solid line corresponds to the ground- 
state of 132 Sn, and the dotted line corresponds to the first excited state. 




Figure 2. Transition strength function \Ma(E)\ 2 assuming discrete levels for 132 Sn. 
See Figure [T] for an explanation of the graphs. 




Figure 3. Transition strength function \My(E)\ 2 assuming discrete levels 132 Sn. See 
Figured] for an explanation of the graphs. 




(Am) 3 

Figure 4. Example of two-particle levels for 132 Sn neutrons (top graph) and protons 
(bottom graph) in the FRDM with LN pairing (solid line). The dotted line is that of 
the Nilsson model with no deformation. 
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(1,1,3/2) 




(n,U) 



(1,0,1/2) 

Levels 
Spherical Shell Model 



Single Particle Levels 



Figure 5. Sample assignment of quantum numbers based on level proximity to levels 
of the spherical shell model [9j. For level e, a contribution from two spherical shell 
model levels is made, and the integration is over the shaded regions in each triangular 
region corresponding to the level. Note that the level height corresponds to the level 
degeneracy, and the widths scale with standard level spacing. 
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Figure 6. Error in calculated /3-decay rates as a function of the known half- life using 
the FRDM [27]. 
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Figure 7. /3-Decay Q-values calculated with the FRDM compared to experimentally 
known Q-values. 




Figure 8. /3-Decay rates of four nuclei as a function of excited state energy of 
individual single-particle excitations. 
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Figure 9. GT transition strength function for various particle-hole states in 
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Figure 10. First-forbidden transition strength function for various particle-hole states 
in 150 Ba. 




Figure 11. /3-Decay rates of the nuclei in Figure [5] as a function of temperature. The 
decay rate of 132 Sn is also plotted as a function of temperature. 




Figure 12. Calculated effective /3-decay rates of r-process nuclei used as a function of 
their mass. The nuclei represented in this figure fall along a line of neutron separation 
energy S„=2.5 McV. The arrows indicate neutron closed shells at N=82 and 126. 




Figure 13. Relative decay rates for various nuclei along the A=162 isobar as a function 
of temperature. In general, the nuclei closer to stability will undergo a larger increase 
in decay rate with temperature. 
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N 

Figure 14. Nuclei included in the current nuclear reaction network (Terasawa et al. 
2001). Nuclei range from stable nuclei, (those on the left side of the indicated region), 
to those along the neutron-drip line (on the right side of the region) . 
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Figure 15. Possible nuclear reactions included in the network calculation. Each arrow 
corresponds to a reaction (and its indicated inverse reaction). Not shown are neutrino 
interactions. Weak interactions are indicated by the lighter arrows. 
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Figure 16. Final r-process abundance distributions for models A, B, and C; shown 
in Figures (a), (b), and (c) respectively. The solar r-process abundance distribution 
(scaled to the figures) is given by the dots. These models do not include /3-delayed 
neutron emission, and only include /3-decays of ground-state nuclei. 
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Figure 17. Comparison of r-process freezeout abundance distributions for models 
B and C in plots (a) and (b) respectively. Solid lines and dashed curves display the 
calculated results with (hot model) and without (cold model) excited-state /3-decays 
respectively. The solar distribution is also shown by the dots. 




Figure 18. Comparison of r-process freezeout abundance distributions for models F 
and G in plots (a) and (b), respectively. Solid lines and dashed curves display the 
calculated results with (hot model) and without (cold model) excited-state /3-decays 
respectively. The solar distribution is also shown by the dots. 




Figure 19. Freezeout r- process abundance distributions for models D and C, as well 
as models E and B in plots (a) and (b) respectively. Both models include excited-state 
/3-decays, while models D and E included /3-delayed neutron emission. 
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Figure 20. Double and single delayed neutron emission probabilities shown in Figures 
(a) and (b) respectively. It can be seen that the nuclei just above the A=130 mass 
region have higher probabilities for two-neutron emission following ,3-decay. 



